Orbital Order, Stripe Phases and Mott Transition in a Planar Model for Manganites 
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Understanding orbital ordered (00) Mott insulating states lies at the heart of a consistent reso- 
lution of the colossal magneto-resistance (CMR) observed in manganites, where its melting induces 
a low-T insulator-metal transition for 0.25 < x < 0.45. Motivated thereby, we study the OO states 
in a planar model for bilayer manganites using DMFT and finite-size diagonalisation methods. We 
derive the correct OO ground states observed in manganites for x = 0, |, |, 4 in exact agreement 
with observations, including the charge-orbital-magnetic ordered stripe phases for x > ^. These OO 



states are exactly shown to be associated with an "alloy" 
on each Mn 3+ site. 

PACS numbers: 71.28+d,71.30+h,72.10-d 



ordering of the d 3 



/d 3y 2_ r 2 orbitals 



(N 



CO 



I 

C 

O 

o 



> 

(N 

in 

m 

o 
o 



Colossal magnetoresistance (CMR) materials have re- 
ceived much attention [1, 2], due to their extreme sen- 
sitivity to minute perturbations [IJ. The parent (cubic 
perovskite) materials are Mott-Hubbard insulators with 
G-type (AF) orbital order of ds x 2_ r 2 /d^ y 2_ r 2 orbitals and 
A- type AF spin order Upon hole doping, x (diva- 
lent ion substitution) in Lai-^Ca^MnOs, for example, 
they evolve through ferromagnetic, orbital ordered (OO) 
Mott insulators with unusual properties 0], to a ferro- 
magnetic metal (FM) at low-T. A transition to a para- 
magnetic insulator (PI), dependent upon cation-dopant 
type, is seen for T > T c [3J. A small magnetic field 
suppresses this FM transition, leading to CMR. These 
phenomena are also seen in bilayer manganites. Fur- 
ther, more "strange" OO states are found in "overdoped" 
(with Ca) manganites. The half-doped manganites show 
a charge (C), orbital and AF order that is very sensitive 
to small perturbations Q (H ext = 5 — 7 T gives a ferro- 
magnetic metal with no CO/OO). The "overdoped" man- 
ganites with x = i §, f , | show extremely stable pairs 
of Mn 3+ C>6 Jahn- Teller distorted stripes having periods 
between 2 — 5a (a=unit cell length); for other values of 
x, a mixture of the two adjacent commensurate configu- 
rations is found Q . For x = 1, CaMn03 is again an AF 
(S = §,i2g) Mott insulator. Finally, the correlated na- 
ture of manganites is shown by dynamical spectral weight 
transfer (SWT) over large energy scales O(4.0 eV) in var- 
ious 0, S, studies as a function of x, T, B ex t-this can 
only result from strong electronic correlations. The im- 
portance of the Jahn- Teller ( JT) coupling is evidenced by 
the large isotope effects .10] and by FM transitions driven 
by O 18 — ► O 16 isotope substitution (see, however, 
Ref. [l2(], where the JT coupling is argued to be much 
weaker than in [l(| 11 1). Thus, understanding CMR is 
inextricably linked to understanding how these strongly 
coupled orbital-spin-charge correlations are modified by 
small perturbations as a function of x. A unified de- 
scription of these unusual observations in one picture is 
a formidable challenge for theory. 



The CMR problem has been extensively tackled in lit- 
erature [i~3L [3] using a variety of numerical and ana- 
lytic (QMC and D = 00) methods, for double exchange 
(DE) models, with/without Jahn- Teller phonons, as well 
as with strong multi-orbital Coulomb interactions with 
static/dynamic JT phonons [l3 ]. For OO states, the full 
multi-orbital Hubbard model has been studied by map- 
ping it to a Kugel-Khomskii (KK) model 15j. However, 
a controlled treatment (semiclassical analysis indi- 
cates an order-by-disorder mechanism) is hard: even the 
type of order is unclear there, and the results sensitively 
depend on the approximations used [l7j ]. 

Here, we take the first step to study the OO, Mott 
insulating phases observed in CMR manganites within 
a 2D, multi-orbital Hubbard model incorporating the 
above-mentioned strongly coupled correlations. Our con- 
clusions apply, with small additional modifications (to be 
treated separately) to bilayer manganites. We show that 
a 2D model suffices to capture the correct OO states 
observed as a function of doping, x, and leave the full 
3D problem for a separate work. Going beyond previous 
studies 14|, [17J , we show how incorporation of the realis- 
tic structure of a single MnO^ layer explicitly in the one- 
electron hopping integrals introduces new, unanticipated 
features, making a qualitative difference to the physical 
results for all x. Further, we show how the "strange" 
stripe-ordered phases in the global phase diagram are 
naturally rationalised from our effective model. 

We start with a model that explicitly includes orbital 
degeneracy of the e g orbitals in manganites 0] , 



H = - t#(at r b jtr + h.c) + u n m n m 

<ij>a,b i,(3—a,b 
iaa' iaa' 

+ H. ]T , (1) 

where the a and b are fermion annihilation opera- 
tors in the doubly degenerate e g orbitals, tfj (a,b — 
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d 3x 2_ r 2, d 3y 2_ r 2) is a 2x2 matrix in orbital space in- 
corporating realistic features of the basic Mn — O per- 
ovskite structure [la ]. U,U' are the on-site, intra- and 
inter-orbital Hubbard interactions, and Jh is the (strong) 
Hund's rule coupling giving rise to the FM state as in the 
usual DE model. Polaronic effects are described by Hjt 
(see below). 

At strong coupling, setting U, Jh >> t 
gives the following effective Hamiltonian, Hq = 



E 



j-ab 



it r Yij(S)(a]bj + h.c) with fi = x,y. Here, 



tf = |[3,V3,V3,1] and tf = |[3, -^3, -y/3, 1] define 
the one-electron hopping matrix for a single manganite 
layer. We now turn on U'. One is effectively dealing 
with spinless fermions, but now with an orbital index. 
Clearly, this model (U' = 0) cannot access the interplay 
between magnetism and 00 in manganites. With U' 
and the JT coupling terms, H becomes 



n ia riib + Hjt 



(2) 



where 7ij(S) is the usual DE projection factor [3j. 

Transform to new variables, c a f = (o+(— l) a \/3b)/^/2, 
c al = ((-l) a y/Ea - b)/V2 with (-1) Q = +1 (a||a;) and 
= —1 (a||y). The c aa transform exactly like d^ x 2_ r 2{] 
), dz y 2_ r 2 (|). This exactly yields a Falicov-Kimball model 
(FKM) where only the c a \ hop; the c Q j_ are strictly im- 
mobile as long as no JT distortions are included. Thus, 

He ff = ~ i 7 l j(S)(c 4 t QT c jaT +/i.c) + J7'^n mT n Mi 

<ij>,a i,a 

+ Hjt = Hp km + Hjt , (3) 

reflecting the correlation between the magnetic and or- 
bital degrees of freedom described above. 

In orbital space, the JT coupling corresponds to ad- 
dition of external fields |18(, Hjt = Q2 J2i( n ia ~ n ib) + 
Q3 J2i( a \bi + h.c). In the rotated basis, this is, 



H 



JT 



Hal 



Hal } 



(4) 

where Q++ = ((-1) Q V3Q 2 - Q 3 )/2 and Q+_ = (Q 2 + 
(— l) a v3<33)/2 are staggered JT distortions which follow 
the orbital (electronic) variables. So H e ff = Hfkm + 
Hjt is a FKM with a local, staggered hybridisation be- 
tween the c a |, c a i a ^ each site. Inclusion of finite phonon 
frequency (MQ 2 {Q\ + Q|)/2) and intersite phonon cou- 
pling terms is required in a full analysis: we have not 
done this here. 

For a half- filled band of spinless fermions, the exact 
solution of H in D — 2 implies an anti-ferro orbital or- 
der of ds x 2_ r 2,ds y 2_ r 2, exactly as required [l9l |. Such a 




FIG. 1: Phase diagram for the transformed multi-orbital 
model with Q++ = at half-filling. The charge Mott (CMI) 
and the band (BI) insulators, both with anti-ferro-orbital or- 
der (AF-OO) are separated by an incoherent, pseudogapped 
metal (M) phase. 



but 



c Q | = d x 2_ y 2 : c a i = d 3z 2_ r 2 there. This would lead 
to an AFOO of d x 2_ y 2 /d 3z 2_ r 2 , at variance with observa- 
tions. Here, such a FKM follows exactly from the realistic 
hopping structure. Moreover, the AF-00 (Mott insulat- 
ing, see below) state is driven by large U' , in contrast 
with band-based scenarios. We note that Yamasaki et 
al. [13] have derived an AF-00 Mott insulator for cubic 
LaMn0 3 (with x = 0) using LDA+DMFT. Our work is 
thus complementary to theirs for x — 0, but goes much 
further, permitting us to study the "exotic" OO states 
for x > h as well (see below). Moreover, given our effec- 



tive FK mapping [2l|, the 00 state(s) are readily under- 



l 3y 2 



FKM has been employed earlier [1J, ll7[ for mangnanites 



stood in terms of an "alloy" ordering of ^33,2. 
orbitals at each Mn site. 

We now solve H e ff = Hfkm + Hjt in d = 00. As 
shown earlier [lit ]. DMFT works surprisingly well for 
the 2D FKM. The FKM with/without Q + _ has an al- 
most exact solution in D = 00 22]. The formalism is 
essentially the same as that used previously, and gives 
very good agreement with QMC results for the same 
model [23| . Keeping U'/t fixed and large, phase transi- 
tions from the Mott insulator with AF-00 to correlated 
(incoherent) metal with no 00, to a correlation-asssisted 
band insulator, again with AF-OO, occur: this is indeed 
borne out in the D = 00 solution, as shown in Fig. [1] 

Given that U' is much larger than in H above, 

we conclude that manganites fall into the CMI class with 
AF-OO, and that the JT terms lead to additional stabi- 
lization of both. Finally, DMFT gives the full, correlated 
spectral functions of the model for arbitrary parameter 
values and band-fillings, at a very modest numerical cost. 
This allows us to study the filling driven Mott transition 
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FIG. 2: (Color online) Partial orbital-resolved (top panels) 
and the total many-body DOS for the AF-OO phase of H e fj 
(see text) within DMFT for U' = 2.6 eV, and staggered 
Q++ = 0.3 eV, Q+- = 0.4 eV, for various band-fillings. Off- 
diagonal components of the spectral fuction are not shown. 
For (n) — 0.9, 0.8, the low-energy pseudogap at u = is 
clearly seen in the DOS. 



from an AF-00 Mott insulator to an incoherent metal 
(see below). 

The relevant DMFT equations were derived ear- 
lier so we do not repeat them here. Since the JT 
terms are staggered, but bilinear in the e g basis, they 
are easily incorporated into the earlier DMFT structure. 
The Green function is now a (2 x 2) matrix in orbital 
space. The staggered, JT "external field" terms im- 
ply an averaging over their orientations, which is car- 
ried out within the DMFT equations to yield the DOS. 
We choose U' = 2.6 eV, Q++ = 0.3 eV, Q+_ = 0.4 
eV as model parameters along with a non-interacting 
DOS for the 2D square lattice with bandwidth, W = 
2.0 eV and variable band-filling, n = (1 — x), in the 
DMFT solution. For n — 1, (see Fig. rjj we obtain 
an AFOO Mott insulator. This is obtained from the 
computed value of D\ a = (— ^) a {{c[ a ^Ciai + h.c)) — 

C(w'<2++.+-) = °- 07 and D *<* = (-!) a ((^aT " 
niai)) = C(w> Q++.+-) = °- 05 ( not shown), ob- 
tained directly from D^ a = — —J (jImG a(T (uj)duj and 
Di a = — —J ImG a ii(uj)duj from the DMFT equations. 
Away from n = 1, the DMFT equations have to 
be supplemented with the Friedel-Luttinger sum rule, 
(n) = J^J2 a ,a ImG aa(^)duj. This is computed 
self-consistently within the DMFT. 

For (n) = 0.9, 0.8, we obtain an incoherent, pseu- 
dogapped, metallic state (see Fig. [2]) with a sharp re- 
duction of local anti-ferro orbital (AFO) correlations 
(Di a — 0.009). Thus, appearance of the doping-driven 



FIG. 3: Two different charge-orbital ordered (COO) ground 
states of the effective Falicov-Kimball model for x — i (left) 
and x — | (right). These exactly correspond to the COO 
states observed in manganites for these x values 0, |f|. 



(FM) metallic state is intimately linked to the "melt- 
ing" of local anti-ferro orbital correlations of the Mott 
insulator with x. The non-FL character of the FM con- 
trasts with what is expected in the FKM with uniform 
hybridisation (V = in the usual FKM with hy- 
bridization), where a correlated FL metal is obtained 
whenever V is relevant (24j. In our model, the staggered 
"fields" produce a low-energy pseudogap, sup- 
pressing FL coherence. Chemical disorder will further 
reinforce incoherence [26| . Given the d-wave character 

of the staggered JT terms (note that both Q++^ have 

components that change sign under a 7r/2 rotation in xy 
plane), as well as the (more important) fact that d-wave 
ground states are obtained near half-filling in a Hubbard- 
like (FKM) model [25[, we predict that this incoherent 
FM- metal phase will exhibit a d-wave pseudogap. 

In contrast to earlier FKM work [l~7l ]. however, 
the ordered, insulating phases in un(doped) manganites 
arise naturally from our model. The checkerboard order 
of d^ x 2_ r 2,d 3y 2_ r 2 corresponds to an AF-00 insulator. 
The exotic bi-stripe states too are naturally predicted 
from the analysis of our FKM. In the insulating phases, 
the "hybridisation" (Qh ) is irrelevant, and the result- 
ing FKM rigorously undergoes phase separation into 
hole-rich (orbital disordered) and hole-poor (orbit ally or- 
dered) phases, as shown by Freericks et al. [2l| by mini- 
mizing the total energy for various x. We have repeated 
their analysis for various x > 0.5. For x = 5, |, |, |, we 
obtain stripe phases with periods 2,3,4,5, as observed 
by Mori et al. (f| using electron diffraction. In Fig. [3l 
we show only the OO ground states for x = ^, |; these 
exactly correspond to those observed in manganites for 
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these hole dopings. 

Given that Mn 3+,4+ correspond to one/zero e g elec- 
tron on each Mn site, the 2D model automatically has 
charge order (CO) of the correct types for these val- 
ues of x. Also, the stripe 00 of pairs of Mti 3+ Oq 
(distorted) octahedra automatically corresponds to a bi- 
stripe charge-order (CO) of e g electrons with the peri- 
odicity determined by x Q. Given the bi-stripe 00 
states, Goodenough-Kanamori-Anderson rules directly 
imply that intersite interactions between the "core" ti g 
spins (S = 3/2) will lead to AF-coupled ladders (Mn 3+ ) 
separated by strips of JT-undistorted (Mn 4+ ) regions. 
Given suppression of e g hopping in an AF "background" , 
these stripe states will be insulators, as observed [3, 
These states will be further stabilised upon inclusion of 
JT terms and longer range elastic interactions. 

This fully corresponds to observations in bilayer man- 
ganites for x > 0.5 [2|. Thus, stripe states in "overdoped" 
CMR exactly result from an "alloy" ordering of a binary 
alloy of Mn 3 +{S = 2, d 4 ) and Mn 4+ (S = §, d 3 ) orbitals 
with d 3x 2_ r 2 , G? 3y 2_ r 2 symmetry. Phase separation/stripe 
phases have long been studied using the FKM (binary al- 
loy disorder model) in alloy physics [27] . Here, we show 
how these phenomena in manganites arise from strong, 
multi-orbital, electronic correlations, which are now ex- 
actly representable as a binary alloy model. Since 00 
states spontaneously break discrete, Ising symmetries of 
H (Eq. {T])), the link to alloy ordering (described within 
an Ising model framework [271 ]) is readily apparent. Re- 
cently [28], 00 phases in a 3D model were derived within 
a static Hartree-Fock approximation. In future, we will 
make contact with these results. 

To conclude, we have shown how consideration of the 
actual multi-orbital structure of the hopping matrix in 
the e g sector within a multi-orbital correlated model re- 
sults in an understanding of the various 00 insulating 
phases observed in CMR manganites, especially in bi- 
layer function of x. These are now under- 
stood simply as an "alloy" ordering of d^ x 2 _ r 2 , d^ y 2 _ r 2 
orbitals, driven predominantly by the inter-orbital cor- 
relations (U 1 ). Our study shows that 00 in overdoped 
(x > 0.5) manganites need not imply very strong JT cou- 
pling, in agreement with [12]: by itself, U' leads directly 
to such phases as a function of x. A moderate JT distor- 
tion will further stabilise these ordered phases. Within 
multi-orbital DMFT, we have shown how an AFOO/F 
Mott insulator turns into a correlated, incoherent, fer- 
romagnetic "bad metal" upon hole doping. This goes 
hand-in-hand with a drop in local AFO correlations. 
These results are fully consistent with indications from 
a host of experiments probing various phases of doped 
bilayer manganites. Interestingly, planar nickelates are 
also modelled by a similar Hamiltonian, and our work 
also naturally explains the OO/stripe phases observed 
there J29j. We expect our analysis to be broadly appli- 
cable to a variety of TMO systems showing a variety of 



OO/magnetic ground states as a function of a suitable 
"tuning parameter" . 
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